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Abstract 

We have been working in many aspects of the problem of analyzing, under- 
standing and solving ordinary differential equations (first and second order). As 
we have extensively mentioned, while working in the Darboux type methods, the 
most costly step of our methods and algorithms of solution is the determina- 
tion of Darboux polynomials for the associated differential operators. Here, we 
are going to present some algorithms to greatly reduce the time expenditure in 
determining these needed Darboux polynomials. Some of them are based on a 
detailed analysis of the general structure of second order differential equations 
regarding the associated differential invariants. In order to perform this analysis, 
we produce a theorem concerning the general form for the differential invariants 
in terms of the Darboux polynomials. 
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1 Introduction 



The differential equations (DEs) are the most widespread way to formulate the 
evolution of any given system in many scientific areas. Therefore, for the last 
three centuries, much effort has been made in trying to solve them. 

In the first approaches to solving DEs, the methods consisted on, for a given 
DE or class of them, trying to use a specific method to deal with that particular 
case. With the cumulation of such experiences, a classificatory system was pro- 
duced where a specific set of rules was designed for certain classes of DEs. By 
the mid XIX th century, this situation had established itself as the way to deal 
with DEs. 

Latter on, more comprehensive methods (non-classificatory) have been de- 
veloped. Just to mention a few, in the second half of the XIX th century S. Lie 
[1, 2, 3, 4, 5], by the development of his continuous groups of symmetries, man- 
aged to unify this whole lot of particular approaches to deal with ordinary DEs 
(ODEs) into a single theoretical framework. Another type of approach we would 
like to highlight is the Darboux type approach [6, 7, 8, 9, 10, 11]. 

In the Darboux type approach front, it is worth mention that a (semi) algo- 
rithmic approach applicable to solving first order ordinary differential equations 
was made by M. Prelle and M. Singer [6]. The attractiveness of the Prelle- 
Singer (PS) method lies in that, if the given first order ODE has a solution in 
terms of elementary functions, the method guarantees that this solution will 
be found (though, in principle it can admittedly take an infinite amount of 
time to do so). The original PS method was built around a system of two 
autonomous first order ODEs of the form x = N(x, y), y = M(x, y) with M and 
N polynomials in C[x, y] or, equivalently, a rational first order ODE of the form 2 
z = dy/dx = M(x,y)/N(x,y). 

We have been working on analyzing and solving systems of first and second 
order differential equations (lODEs and 20DEs, respectively) from a numerical 
point of view [12], using Lie methods [13, 14] and Darboux type approaches 
[15, 16, 17, 18, 19, 20, 21, 22, 23]. For this latter class of methods, we have been 
developing (semi) algorithms to deal with classes of ODEs. In these algorithms, 
one fact has been always present: the most (computationally) costly step is the 
determination of the associated Darboux polynomials. 

Based on that realization, here we will be focused on speeding the process of 
finding Darboux polynomials for a class of ODEs of our interest. In particular, 
in this paper, we will talk about a class of rational 20DEs. 

This finding of the Darboux polynomials (quicker than before, via other meth- 
ods) is useful on its own right. But, in particular, its very compelling to the class 
of Darboux type methods we have been working with. To exemplify this, we will 
use the algorithm developed by us in [22]. 

The paper can be summarized as follows: In the next section, we will briefly 
introduce the main features of the method we will use as model for the methods 
benefitted from our here presented method to speed the finding of the Darboux 
polynomials. In section (3), we will present a part of the method where one can 
find, in a very straightforward manner, the Darboux polynomials for a class of 
20DEs by inspection. Next, in section (4), we will present the part where a 
deeper analysis of the structure of the 20DE is needed to extract the Darboux 
polynomials. In each of the above mentioned two sections, we will present ex- 
amples of the application of our proposed method. Finally, we will present our 



2 From now on, we will use that dy/dx = z. 



conclusions. 



2 Summary of the Darboux type approach 
we will use as model for a the benefited ap- 
proaches 

In [15] we developed an extension of the Prelle-Singer method [6] and, in that 
paper we proposed to use an unknown function (that we called S) in order 
to make the 1-form 3 <f)(x, y, z) dx — dz proportional to an exact 1-form. In 
[21] we constructed a semi-algorithm to determine the 5-function for an 20DE 
presenting an elementary first integral. In [22], we built on a new, more efficient 
algorithm, and established a set of theoretical results that makes the basis of 
the algorithm. The work here presented will, in the majority, be concerned in 
improving one step of the algorithm: namely the finding of Darboux polynomials 
(that is an essential part of the procedure). 

In [22], we use a variation of the idea of the function S. Basically, we use (in 
order to span the space of null 1-forms) that: 

If the 20DE is given by: 

d 2 y x M(x,y,z) . . 

where M and N are polynomial functions of (x, y, z). 
We then construct a differential operator D, 

D = N d x + z Nd y + M d z , (2) 

extracted from the second order ODE, and their corresponding Darboux poly- 
nomials and co-factors will be the building blocks of the integrating factor and, 
ultimately, of the differential invariant for the 20DE. 
We can write (1) in the form 

dx dz 

~N = M ^ 

and since z dx = dy, we have that the 1-forms a and j3 defined by a = M dx— N dz 
and j3 = zdx — dy are null over the solutions of the ODE (1), i.e., over the 
solutions, 

a = M dx - N dz = 0, ft = zdx - dy = 0. (4) 
From these results we have the following: 

If I(x,y,z) is a first integral of the ODE (1) then the 1-form dl is a vector in 
the subspace spanned by the 1-forms a and j3 defined above, i.e., 

dl = ra + s/3 (5) 

where r and s are functions of (x, y, z). 

From (5) we have 

I x dx + I y dy + I z dz = r (M dx — N dz) + s (zdx — dy) (6) 



3 This 1-form is associated with the rational 20DE z' — <p(x,y,z), where is a rational functions 
of (x,y,z). 



implying that 



I x = rM + sz 



h = -rN 



(7) 



Therefore, if we determine, r and s we can find I via quadratures using the 
equation: 



I(x,y,z) 
d 



Jl x &X + J[l y --l y JlA.r 



I x dx 



dy 



dy + 

dz 



(8) 



In [22], in order to generate an operational semi-algorithm, i.e. being able 
to determine I x , I y and I z , we have used the class of ODEs where s/r (r and s 
defined by equation (5)) is a rational function of (x, y, z) ( i.e., s/r = P/Q where 
P and Q are polynomials that do not have any common factors). 

Also in [22], we have demonstrated many results, including two theorems, 
that allow us to construct a semi-algorithm to find elementary first integrals 
of a class of rational second order ODEs (via a Darboux-type procedure), thus 
generalizing the correspondent method developed by Prelle and Singer for first 
order ODEs. Brutally summarizing the important results for us here, we have 
found that: 



We can write: 

^/Q=n-r^^=E-*^ M =E™^ o) 

r- It — Vi 

where R is the integrating factor and V{ are irreducible Darboux polynomials (in 
(x,y, z)) of the D operator and the gi are the corresponding co-factors. Then, 
from the compatibility conditions, (I xy = I yx ,I xz = I zx and I yz = I zy ), we can 
write: 

P = p ^2 m i9i = ~D[P] - Q(N My - MN y ). (10) 



and 

Q (°W 



u ) =Qj2 m i9i = -D[Q]-P-Q(N x + N y z + M z ) (11) 
Ultimately, these ideas led to: 

dI = R [(MQ + zP)dx + {-P)dy + {-NQ)dz] , (12) 

and 

I x = R(MQ + zP) 

^ = -RP 

I z = -RNQ (13) 

In [22], equations (10,11) were the basis of our procedure. So, by determining 
the Darboux polynomials and solving them, we would have found P, Q and R, 
thus enabling us to use equations (8, 13) and find the invariant. 



The great advantage of our method is that it converts the search for first 
integrals into solving (essentially) first degree algebraic equations (in the same 
way as Prelle and Singer did for first order ODEs). As a consequence of this, 
our approach is capable of analyzing the integrability regions for the seconde 
order ODE (for the case where it presents undetermined parameters). To do this 
we have to 'add' some (or all) parameters as variables into the algebraic system 
we have to solve. The solutions exists for values where there is an integration 
possible 4 . Furthermore, the procedure is semi-algorithmic and, given enough 
time, if the solution exists, it will find it (as does the Prelle-Singer approach for 
first order ODEs). 

It is worth to point out that another positive aspect (that will become clear 
just a few steps ahead, on the next section) is that we use (in our method 
presented in [22] the differential operator given by equation (2), instead of the 
(probably) more commonly used one: 

D = d x + zd y + (/)d z , (14) 

where (p has the same meaning as in equation (1). 

The Darboux polynomials for the operator given by equation (2) are used (as 
explained above) in the algorithm (together with their associated co-factor) in 
order for us to find the differential invariants for the 20DE in question. Next, 
let us try to produce ways to find them more quickly. 

3 Darboux polynomials as straightforward 
factors in the numerator and-or in the denom- 
inator 

In this section, we will extract information regarding the Darboux polynomials, 
correspondent to the operator (2) related to the 20DE being studied in a very 
direct way. This, although a simple procedure, will prove essential to solve (or 
at least reduce) some ODEs. 

As mentioned in the previous section, the fact that our method [22] uses the 
differential operator given by equation ( 2 ) is very advantageous. One can see, 
just by inspecting this operator, that if the denominator of the 20DE (1), i.e. N, 
is a function of (x, y) only, it will be certainly a Darboux polynomial of the D- 
operator defined in (2). Actually, for that matter, any factor of N that depends 
only on (x, y) will do. Analogously, for the case of the numerator in (1) but, this 
time, regarding functions of {z) only. The numerator or any factor of it (that is 
a function of z alone) is a Darboux polynomial of the D-operator (2). 

In the next two sub-sections, we will study each of these cases. 

3.1 Darboux polynomials as functions of (x,y) only 

In order to analyze this case, let us first re-write the D-operator (eq. 2): 

D = Nd x + zNd y + Md z , 
In [19, 21] we had already introduced this procedure 



As explained above, the case of interest here is the one where N is a function 
only of (x, y) or it has factors that are so. To cover all the cases, consider the 
following general form for N: 



k 

N=(l[n i (x,y))M(x,y,z) (15) 
i=i 

where rii(x,y) and ftf(x,y,z) are polynomials in (x,y) and (x,y,z) respectively. 
It is easy to see that this general form for N covers all the cases of interest: 
namely, if the there is no factor that is a function of (x, y) one can consider that 
the product (Yli ni(x,y)) = 1 and there is only M(x,y,z) left. This case will 
not be of interest here. On the other hand, if the whole iV is a function of (x, y) 
only, this can be translated by saying that the product (Ylf rii(x , y)) = n\(x,y) 
and J\f(x,y,z) = 1. Of course, the case where there are a few factors is covered 
with k > 1. 

These cases are listed bellow with the general expressions for the Darboux 
polynomials and the associated co- factors: 



case 1: where the whole N is a polynomial in (x,y) 

N = m(x, y) -> D[m] = m {-^) +n lZ (-^) (16) 

So, if v\ = n\ is a Darboux polynomial of D, the associated co-factor is: 

,dni. ,dni. 
* = ^ + Z ^ (17) 



• case 2: where there are one or more factors in N that are polynomials 
without z 



N = (jJ ni (x,y)^Af(x,y,z) ^ 

D[n 3 ]= = ^f[n t (x,y^Af(x,y,z) ^) + z (18) 

where 1 < j < k. So, if we consider v\ = n,j as a Darboux polynomial of 
D, the associated co-factor is: 



91=^ fl V) j *{x, y, z) ((^) + z (^)) (19) 

3.2 Darboux polynomials as functions of (z) only 

This situation is analogous to the one just dealt with. So we will begin the 
exposition in a similar fashion, 

Again, let us reffer to the D-operator given by eq. 2. 

The focus now is when M is a polynomial only on (z) or it has factors that 
are so. To cover all the cases, consider the following general form for M: 



M = (J[ mi (z))M(x,y,z) (20) 
i=i 

where rrii(z) and Ai(x,y,z) are polynomials in (z) and (x,y,z) respectively. 
Again, one can see that this general form for M covers all the cases of interest: 
namely, if there is no factor that is a function of (z) one can consider that the 
product (nf nT-i(z)) = 1 and there is only M.(x, y, z) left. If M is a polynomial on 
(z) only, this can be translated by saying that the product {Yli rrii(z)) = ml(z) 
and M(x,y,z) = 1. Of course, the case where there is a few factors is covered 
with k > 1. 

These cases are listed bellow with the general expressions for the Darboux 
polynomials and the associated co- factors: 

• case 1: where the whole M is a polynomial in (z) 

M = m 1 (z)^D[mi]=mi(^) (21) 
So, if v\ = mi is a Darboux polynomial of D, the associated co-factor is: 

m = (-£-) (22) 



• case 2: where there are one or more factors in M that are polynomials 
without (x, y) 

M= (jlmi(z^j M(x,y,z) -> D[mj] = (jlm^ M(x,y,z) (^f) 

(23) 

where 1 < j < k. So, if we consider v\ = nij as a Darboux polynomial of 
D, the associated co-factor is: 

91= ( U Tm(z)\M(x,y,z) {^f\ (24) 

3.3 Examples 
3.3.1 first example 

Let us present an example that brings Darboux polynomials from M (deppending 
on {z)) and N (deppending on (x,y)): 

d 2 1/2 (2z + 3)(3zy 2 + z + x-y 3 -y-l) 

dx 2 x — y 3 — y 

Using what we have been learning, we can see that we have two Darboux poly- 
nomials: one from M and another from N. Bellow we will display both with 
their corresponding co-factors: 



v\ = 2 z + 3 — > g\ 
v 2 = x - y 3 - y ->• g 2 



= -2 (3 zy 2 + z + x - y 3 - y - l) 

= 2-6zy 2 -2z (26) 



Using the method briefly described in section (2) one conclude that, for this 
ODE, we have the following results for the parameters and functions needed to 
find the differential invariant for the ODE: 



P = (3y 2 + l)(3 + 2z) 
Q = 1 

R = 1 

(3 + 2z)(x-yZ-y) 

and, using (13) and (8), we finally find: 

I = x-l + \n( 2Z + 3 V (28) 
\x - y 4 - yj 

It worth mention that the presence of the Darboux polynomial (x — y 3 — y), 
with the cubic term, makes the regular process of determinig it very "expen- 
sive" in time expenditure. After applying the method here presented, which 
very quickly determined the needed Darboux polynomials, the algorithm we in- 
troduced in [22] finds the results (27) and (28) almost instantly. 

3.3.2 second example 

This second example shows a case, where one of the Darboux needed to deter- 
mine, via the results presented on section 2, is easily determined using the regular 
set of procedures (and this one is not apparent on the format of the 20DE in the 
fashion we are advocating here) while the other, that would be very demanding 
on time, is determined analyzing the denominator of the 20DE. Let us present 
the ODE 5 : 



(27) 



d 2 y (2 Cbxy + 2Cb- xb) z 2 

dx 2 ~ 2 {xy + 1) {Bax + Bbz + Bey + b) + 

(-yb + 2 Cax 2 y + 2 Cax + xcy + 2c + 2 Ccy 2 x + 2 Ccy - x 2 a) z 
2 {xy + 1) {Bax + Bbz + Bey + b) 
yax + 2 a — y 2 c 
2 {xy + 1) {Bax + Bbz + Bey + b) 

For this differential equation, if we we run the regular procedures in order 
to find the Darboux polynomials, it is quickly found one such polynomial. But 
one needs more information than that in order to be able to build an integrating 
factor for (29). Using the analysis presented on section (2) it is easy to see 
that x y + 1 should be a Darboux polynomial for the corresponding D-operator 
(2). This combination of approaches generates the following list of Darboux 
polynomials and co-factors: 



(29) 



v\ = ax + cy + zb — > gi = (2 zxcy + 2 yax + 2 a + 2 zc) B + 

(-2 bzxy - 2 bz) C + xbz + yb 
v 2 = {xy + l) -»• g 2 = 2 {xy + 1) {zx + y) {Bax + Bbz + Bey + 6)(30) 

5 Incidently, althought it is not the point here to compare every method of solving ODEs, this 
particular ODE is not solved by the power methods and algorithms implemented on the powerful 
Maple solve (release 10) 



Using this, one can determine R, P and Q corresponding to (29) (see section 
(2)): 



P = (2ax 2 y + 2ax + 2bzxy + 2bz + 2cy 2 x + 2cy}C-x 2 a-xbz + 

+xcy + 2 c 
Q = 1 

R = "7 7T7 — : 7. (31) 

(xy + 1) (ax + bz + cy) 

and, using (13) and (8), we finally find: 

^ ^ , / ax + bz + cy \ . . 

I = Bz + Cy + In (—j===^-) ■ (32) 

Analyzing the expression for R, we see that we needed two Darboux poly- 
nomials in order to "build" it: namely v\ = ax + cy + zb and i>2 = (xy + 1 
leading to R = J2i V T % R = — (xy + 1) _1 (ax + bz + cy)" 1 . We can also ob- 
serve that, even in this simple case where one of the Darboux is of order one 
and the other of order two, the solution of the ODE can elude many powerful 
techniques (see footnote) . The finding of the "second" Darboux polynomial via 
the technique here suggested proves to be essential to render the Darboux type 
approach presented in [22] practical. 



4 Darboux polynomials from a deeper anal- 
ysis of the 20DE 

In this section, we will deepen our method through a more detailed analysis of 
the structure of the differential equation we want to solve (or, at least, reduce). 
Let us do that beginning with the following: 



Consider that 



I = I(x,y,z) 

is a differential invariant corresponding to a certain 20DE 

d 2 y M(x,y,z) 

-7-2 = <P\ x i y> z ) = T77 7 

dx A N(x,y,z) 

So, by using equation (2), one can write: 



D[I] =04 N 8 X I + z Nd y I + M d z I = ^ 



N I z y > 



where I a = d a I = 91 , a = x, y or z. 



We may ask in what sense this is going to help our quest for Darboux poly- 
nomials. Actually, we now have a great amount of information regarding the 
general structure of the differential invariant that, as we shall show now, will 
help us. 



From [6], we know that the differential invariant is of the form: 

I = W + hi(Wi) (34) 

where Wq and W\ are algebraic functions of (x,y,z). 

Since we are interested in rational 20DEs, one can easily see (from equations 
(33) and (34)) that we would have a relation between Wq and W\. So, from our 
experience in dealing with that question, we do not compromise the generality 
of the method much further if we use the differential invariant in the form: 



A(x,y,z) 
V(x,y,z) 



+ In 



B(x,y,z) 
C{x,y,z) 



(35) 



where A (x, y, z) ,V (x, y,z) , B (x, y, z) and C (x, y, z) are all polynomials. 

Further still, as we will demonstrate below, if equation (35) applies , we 
have that T> (x, y,z) ,B (x, y, z) and C (x, y, z) are Darboux polynomials of the 
D-operator (2). This knowledge will prove essential in the producing of new 
possibilities for determining the Darboux polynomials for a given ODE. 

4.1 Darboux polynomials in the build-up of differen- 
tial invariants 

Actually, in order to demonstrate the above mentioned results, let us redefine a 
different differential invariant such that: 



e 7 = e 



B 



(36) 



\C(x,y,z)J \C, 

where we have dropped the explicit notation of the (x, y, z) dependence but it is 
still there. 

Since X is an invariant, we have that: 



D[l] 



= 



(37) 



So, by using eq. (36), one gets: 



D 



+ 



(?) 



D 



(g)e(*) 



= 



that finally results on 



D 



(?) 



+ D 



(38) 



(39) 



Considering that A,B,C, T> are all polynomials and that we can therefore 
write them in terms of irreducible polynomials p«, we can use the following: 



(f)-n# 



(40) 



in turn, that leads to: 



D 



UkPk 



c k 



E 



'cjD\pj] \ Ej CjI>[Pj] (lI/,, ; />/, 



Pi 



Ilipl 



Using equations (41) and (39), we get: 



(41) 



Ili Pi 



i Wi / \ * / 

Equation (42) can finally be written as: 



= 



(42) 



E c ^N IIp* + \Ii pi ) D 
E^NIIlpfc] + (llPi) D 



Ki + 
Since ifi is a polynomial, so is K2 => 



#2 



(43) 



K 2 



lift 



.d[.A]p-.A£>[p] 

P2 



p(x,y,z) (44) 



where p(x, y, z) is a polynomial. 

Now, remembering that pi are all irreducible, independent polynomials, we 
have that: 



(JhPi 
V 2 



^ polynomial 



(45) 



That situation leads to two possibilities: 

4.1.1 (IliPi) does not have common factors with V 2 

In this scenario, from equations (44 and 45), using equation (39) we have: 



D[A]V-AD[V] 



V 2 



= D 



is a polynomial 



D 



(?) 



-D 



Pi(x,y,z) 



(46) 



where p±(x,y,z) is a polynomial. 



Let us analyze these results: 

Consider the second (the one to the right) equality of equation (46), 



D 



~A 
V 

D[A] _ AD[V] 
V V 2 
AD[V] 
V 



= -pi(x,y,z)=> 
= ~Pi(x,y,z)=> 

- -p 1 (x,y,z)V-D[A] (47) 



Since the right-hand side of (47) is obviously a polynomial and, by construction, 
^ can not be simplified (otherwise it would have been already), we may conclude 
that: 

- is polynomial =4> T> is a Darboux polynomial (48) 

Actually, we can conclude more than that. If we write V as Y\^ q d^ J , where dj 
are irreducible polynomials and kj integers, we have 

D[V]_D[UZi^} _P* U D[<k] 



E**^- ( 49 ) 



=i 



If we multiply (49) by ITj=2 d i > we S e ^ 

n ^ = h (ft ^ + X> ( . ft/;) m- (so) 

Since he left hand side of (50) and the second term on the right hand side 
of (50) are polynomials, we may conclude that k\ (l~[j=2 d j) D[d\]/d\ is also a 
polynomial. Considering that the d's are independent (by construction), the 
product Y[j=2 d j can n °t cancel d±. Therefore, we can conclude that <Zi|.D[<2i]. In 
an analogous way, we have that di\D[di],i = 2 ■ ■ ■ n q . 

So, what we can actually conclude is that V is formed by factors that are 
Darboux polynomials! 

For its importance and the fact that this result will be used below, let us 
state it as a lemma: 



lemma: // a polynomial is written as ]Tj 9 <fr/ , where are irreducible poly- 
nomials and kj integers, and = A, where A is a polynomial and D is a 
differential operator, with polynomial coefficients, we have that: 

4>i\D[4>i],i = l---n q 



Let us now concentrate on the other part (the initial one) of equation (46) 
and, after some algebra, one gets: 



D [B] C - D [C] B = p 1 (x,y,z)CB 
D [B\ D [C] 



B C 



= Pi(x,y,z) 



Pl (x,y,z)B-D[B] = -^P-B (51) 

Regarding the last line of equation (51) above, we can see that the left-hand 
side is obviously a polynomial. Since |?, by construction, is not a polynomial 
(otherwise it would heve been simplified already) one can conclude that: 

— -f^- is a polynomial (52) 

Using this information and the third statement of equation (51) it is straight- 
forward to conclude that 

—p is a l so a polynomial (53) 

So, using the above lemma, one may conclude that: If B is written as B = 
n™ 6 bj J and C as C = YYJ C Cj J , where bj and Cj are irreducible polynomials and 
Kj and A integers, we may conclude that: 

b i \D[b i ],i = l---K q (54) 

and 

Ci\D[ci],i = l---\ q (55) 

So, summarizing the analysis above, for the case where (ILPi) does not have 
common factors with V 2 we have that V, B and C are all built from Darboux 
polynomials. 

4.1.2 (YliPi) does have common factors with V 2 

In order to analyze this case, let us start by writing: 

■ M \ / N \ I M 



n«) = fn«) = fn«)-( n * H- c 

V i J \i=l J \i=l J \i=JV+l / 



(56) 



and 



N \ / N 



V = 11^ • lift = 6 -V (57) 



,i=l / \i=l 



where 




(58) 



Now, from equations (44,56,57,58), we can write: 



lift 



c 



D[A]V-AD[V] 
D[A]V-AD[V] 

ve 



P(x,y,z) 



D[A]V-AD[V] 



VOr, 



(59) 



where p(x, y, z) is a polynomial. 

After a little algebra, equation(59) leads to: 

(D[A]-(A^ = p(x,y,z)9 



(A 



V 
V 



polynomial 



(60) 



Since, by construction, A does not have any factors in common with V and, 
analogously, ( is the part of (IliPi) that "is not" on V, one may conclude that: 



DPI 
V 



b(x, y, z) = polynomial 



(61) 



So, again, by using the lemma proved above, If we write T> as 9 dj J , where 
dj are irreducible polynomials and kj integers, we have: 

di\D[di],i = 1 • • • n q 



Therefore, as in the case of the previous subsection, T> can be written as a product 
of Darboux polynomials. 

Let us now use equations (39,40). Using these, one can write: 



Hi Pi 



+ D 



+ D 



^ c 3 D \Pj] , D 



Pj 



= 

= 
= 



(62) 



Now, remembering that, we are considering the case where T> has factors in 
common with (IliPi) (equations (56,57)) and that we have just proven that V 
can be written as a product of Darboux polynomials, we can write: 



N 

E 



c jD[Pj 

Pj 



M 

E 

j=N+l 



c j D [Pj 
Pj 



+ D 



pi (polynomial) 



D 



Y,f=N+lCj D \Pj\ [Uk^jPk 



+ Pi 



Hi^N+l Pi 



(63) 



where the last line was obtained by re-writing the summation with a single de- 
nominator. 

Finally, one can conclude that: 



D 



D 



M 

n k 

i=N+l 



+ Pi 



M 



2^j=N+l c 3 



D[Pj] (Uk^jPk) 



Pi 



M 



= -Pi II Pi] - J2 c i D \Pj\ \{Pk 



(64) 



> i=N+l 



j=N+l 



A direct inspection of the equation above allows for the conclusion that: 

p2 = polynomial 



D 



M 

n Pi 

,i=N+l 



(65) 



Expanding the left-hand side of the equation above and considering that, by 
construction, (ili^jv+i Pi) does not have any common factor with V, we have 



D[A]V-AD[V] 



V 2 



M 

n « 

K i=N+l 



P2 



D 



p3 = polynomial 



(66) 



So, using this, and the reasoning following equation (46), we may conclude that, 
for the situation being studied here (i.e., (IliP«) does have common factors with 
V 2 ), we also have that 



D 



(?) 



®(x,y,z) 



polynomial 



D[B] . , . , D[C] . 

— — is a polynomial and — — is a polynomial 



(67) 



So, by the lemma proved above and by equations (48,67), we may infer that: 

If we write V as JX™ 9 > B as jXj 1 ' 1 ^j 3 an( i ^ as JJj q , where dj , bj and cj 
are irreducible polynomials and kj, Kj and Vj integers we have: 



bi\D[bi],i = l---n q , bi\D[bi], i = 1 ■ ■ ■ m g and bi\D[bi], i = 1 ■ ■ ■ n q (68) 

So, also for the case where (IliPi) does have common factors with V 2 we have 
that V, B and C are built from Darboux polynomials. 

So, summing up the results for sections (4.1.1,4.1.2), we conclude that for a 
differential invariant given by equation (36), the following result is demonstrated: 



Theorem: For a differential invariant of the form X = e(^) (u)' w ^ ere A,T>, 
B andC are all formed by products of independent, irreducible polynomials (i.e., 
the fractions can not be simplified further) such that T> as Y\ r - q d- j , where dj are 
irreducible polynomials and kj integers, B as Yl" 1 " 1 b - , where bj are irreducible 



polynomials and \j integers and C as JT^ 9 Cj > where Cj are irreducible polynomials 
and bj integers, we have that the dj,bj and Cj are all Darboux polynomials of the 
operator D = N d x + z Nd y + M d z 

So, combining all the above discussion, equation (35) leads to: 

(B x + zB y ) CV 2 + (-C x - zC y ) BV 2 + (A x + zA y ) BCV + {-AV X - zAV y ) BC 
® ~ -A Z V CB + AV Z CB - V 2 B Z C + V 2 BC Z 

(69) 



As we shall see in the next sections of the paper, the information contained 
on (69) will allow for finding "hidden" Darboux polynomials. 

4.2 Darboux polynomials as functions of (z), present 
on the denominator 

One might ask what is the point of the title for this subsection. As we have 
shown previously, from the denominator of the ODE (or from its factors), if it 
(they) is (are) function(s) of (x, y) only, we can extract Darboux polynomials 
straightforwardly. We have shown in section (3.1) that this information might 
prove crucial to reducing the ODE in question. So the question that one may 
ask is if it is possible to extract information regarding Darboux polynomials 
depending on z from the denominator of the ODE. Using (69) one can. 
Let us consider the denominator of the ODE as written on (69): 

N = -A Z VCB + AV Z CB -V 2 B Z C + V 2 BC Z (70) 

Of course, that is, in principle, the general expression for the denominator of 
the ODE. Surely, it can happen that, for some particular combination of A,T>,B 
and C, some factor of (70) will be also a factor on the numerator for the ODE 
and some cancellation will occur. But, disregarding this possibility for now, let 
us particularize the analysis a little. 

Let us consider the particular situation where: 

V = k(constant) 

A = A(x,y)^A z = (71) 
For this case, equation (70) becomes: 

N = -V 2 B Z C + V 2 BC z ^k 2 {-B z C + BC Z ) (72) 
Consider the case where: 

B{x,y,z) = B(z) + b(x,y) 

C{x,y,z) = C{z) + c{x,y) (73) 

where B,C,b and c are polynomials. 

So, combining these results (eqs. (71) and (73)), we may finally conclude 
that: 

N = k 2 (-B z C + BC z )^-k 2 ^B(z)^C(z)-k 2 ^B(z)y(x, y ) + 
+k 2 (±C (*)) B (z) + k 2 (±C (*)) b (x, y) (74) 



It seems that the situation has not improved. But, take heart, if we further 
consider that 

B(z) = C{z) = K{z), (75) 

we finally find that: 

N = -k 2 ^B(z)y(x,y) + k 2 ^ z C(z)y(x,y)^ 

-+fc 2 (^/C(z)) (b(x,y)-c(x,y)) (76) 

Actually, a little variation of this idea could broaden the range for the appli- 
cation of it. Consider that, instead of equation (75) we had: 

c B B{z) = c c C{z) = c B c c IC(z), (77) 

where cb and cq are constants. In turn, that will imply that 

B z = c c K-z 
C z = cbK, z — >■ 

= k 2 ^lC(z)yccHx,y)-c B c(x,y)) (78) 

If we use cb = cc = in (77) we recover (75). Thus, using (77) generalizes (75). 
This can be very useful, as we shall see on section (4.4). 

So, the main trust of the approach suggested in this section can be summarize 
as an algorithm in the following steps: 

Steps of the Algorithm 

1. Inspecting the denominator of the 20DE, analyzing equation (76), deter- 
mine candidates for the functions )C(z), B(x, y, z) and C(x, y, z) (eq. (73)). 

2. Considering that the conditions of sub-section 4.1 are met, these functions 
would be then, each, Darboux polynomials. Verify that. 

3. If the above verification proves true, run the algorithm introduced in [22]. 
An example of the application of these ideas will be presented in section (4.4). 

4.3 Darboux polynomials as functions of (x, y), present 
on the numerator 

Here we are going to do something similar to the procedure shown above (where 
we have extracted information enabling us to determine Darboux polynomials 
depending on z from the denominator of the corresponding ODE). We are going 
to study the possibility of extracting Darboux polynomials (deppending on (x, y)) 
from the numerator of the ODE. 

Let us analyze equation (69). Restricting ourselves to the case: 



T> = k (constant) 

A = A{z) -> A x = A y = 0, 



(79) 



the numerator for the ODE (see (69) becomes: 

M = (B x + zB y ) Ck 2 + {-C x - zC y ) Bk 2 (80) 

Surely, it can happen that, for some particular combination of A,V,B and 
C, some factor of (80) will be also a factor on the denominator for the ODE and 
some cancellation will occur. But, disregarding this possibility for now, let us 
particularize the analysis a little further. 

Consider the following: 

B(x,y,z) = P{x,y) + \>(z) 

C{x : y,z) = 7(z, + (81) 

where /3,7, b and \ are polynomials. 

Following the lead from sub-section (4.2), equation (75), considering: 

(3(x,y) = j(x,y) = 9(x,y), 

after a little algebra, from equation (80) one gets: 

M = (B x + zB y )Ck 2 + {-C x - zC y )Bk 2 -»• 

-> k 2 (e x + z6 y )^{z)-\>{z)) (83) 

Again, we may summarize the procedure here introduced to try and find 
Darboux polynomials, as functions of (x,y), from the numerator of the ODE 
under consideration, as an simple algorithm: 

Steps of the Algorithm 

1. Inspecting the numerator of the 20DE, analyzing equation (83), determine 
candidates for the functions 9(x,y), \\{z) and \>(z). 

2. Considering that the conditions of sub-section 4.1 are met, these functions 
would be then, each, Darboux polynomials. Verify that. 

3. If the above verification proves true, run the algorithm introduced in [22]. 

In the next section (4.4), we will present an example where this equation is 
put into use and provides the finding of Darboux polynomials. 

4.4 Examples 

Although, as mentioned above, it is not the point of this paper to compare the 
capabilities of the method presented on [22] with other methods and algorithms 
that deal with ODEs, it is interesting to mention, in the sense of displaying the 
practical usage of the ideas here presented, that both examples on this section 
can not be solved by the powerful methods implemented on the Maple symbolic 
basin (release 10). 



(82) 



4.4.1 third example 



Here we are going to exemplify the usefulness of the ideas and equations devel- 
oped in section (4.2). 

Consider the following ODE: 



d 2 y —cz 10 a — 6 z 6 dy 5 a + (—be + ady 6 — cbx) z 5 — 6 zdy 5 bx + dy 6 bx + bdy 6 
dx 1 ~ (-5z 4 ) (ady 6 + cbx) 

(84) 

Comparing this above equation with equation (76), we can try and deduce two 
Darboux polynomials from z 4 (ady 6 + cbx) . 

The most obvious "guess" is that, for example: 

k = 1 

K(z) = z 5 

b(x,y) = ady 6 

c(x,y) = —cbx (85) 

Unfortunately, this does not lead to actual Darboux polynomials. So it seems 
that the analysis presented on section (4.2) does not help us here. But, as pre- 
viously advertised, if we use equation (77), we would be lead to make a different 
"guess" and that might result successful. Trying 



k 


= 1 


K(z) 


= z 5 


B(z) 


= alC(z) 


C(z) 


= cJC(z) 


b(x,y) 


= dy 6 


c(x,y) 


= —bx 



(86) 

we find a couple of Darbox polynomials. We present them and the corresponding 
co-factors bellow: 

v\ = —cz 5 + dy 6 g\ = 5 z A (ft zdy^a + z 5 ac + cbx + bej 

v 2 = z 5 a + bx g 2 = -5 (-6 zdy 5 a -bc + ady 6 - z 5 acj z A (87) 

As usual, using these results, running the methods and algorithms introduced 
on [22], we can find: 

P = -6dy 5 (z 5 a + bx) 
Q = 1 

R ~ (-cz 5 + dy 6 ) (z 5 a + bx)' ^ 
and, using (13) and (8), we finally find: 

I = -x + \n (-cz 5 + dy 6 ) - In (az 5 + bx) . (89) 

From (88), it is easy to see that we needed the two Darboux polynomials we 
have found using the technique displayed on section (4.2) in the sense that the 
integrating factor R is constructed from both of them. 



4.4.2 fourth example 



In this section, we are going to exemplify the use of the results listed on section 
(4.3). Let us get to work by examining the following ODE: 

d 2 y _ - (az 2 - z -2az + a) (Azy 3 + 3bx 2 ) 

dx 2 ~ (bx 3 + y 4 ) (3 z 4 a - 3 z 3 - 6 z 3 a + 3 z 2 bx 3 + 3 z 2 y 4 + 2 az 2 + z + a - bx 3 - y 4 + 1) 

(90) 

Using the techniques suggested on sections (3.1) and (3.2), we can deduce that 
(bx 3 + y 4 ) and (az 2 — z — 2 az + a) are Darboux polynomials of the D operator 
(equation (2)) corresponding to the ODE (90). Indeed, this is true and very quick. 
The question is: is that enough to find an integrating factor (and, consequently, 
a differential invariant) for (90), using the Darboux method we are using here 
[22] The answer is no. So, let us ask section (4.3) for help. 

Basically, we want to compare equation (83) with the numerator of the ODE 
(90). 

k 2 (6 X + z e y ) (l](z) - \>{z)) (az 2 - z - 2 az + a) (4 zy 3 + 3 fcc 2 ) 

(91) 

Surely, there are a lot of possible sets of \>(z),\\(z),k and 9(x,y) to satisfy 
equation (91). A satisfactory answer (one producing Darboux polynomials in 
enough numbers to find a differential invariant for the ODE under study) is 
found for: 



k = 

\>{z) = 

\\(z) = (az 2 — z — 2 az + aj 

0{x,y) = [bx 3 + y A ) (92) 

With this solution to (91), we find the following two Darboux polynomials 
and corresponding co-factors: 

vi = (bx 3 + y 4 ) -»• 

gi = (3 bx 2 + 4 zy 3 ) (3 z 4 a - 6z 3 a-3z 3 + 3 z 2 bx 3 + 3 z 2 y 4 + 2 az 2 + z + a - bx 3 - y 4 + 1) 
V2 = (bx 3 + y 4 + az 2 — 2 az + a — z) — > 

92 = - (3 bx 2 + 4 zy 3 ) (6x 3 - 3 z 2 bx 3 - 1 + y 4 - 2 a + 2az - 3 z V) 

(93) 

using the method briefly explained in section (1), we calculate 



P = 4y 3 (az 2 - z-2az + a) 
Q = 1 
R 



(bx 3 + y 4 ) (bx 3 + y 4 + az 2 - 2az + a - z)' 



(94) 



and, finally: 



5 Conclusion 



From our experience working on the Darboux type methods to deal with ODEs, 
we concluded that, on all the methods and algorithms, the necessary step of 
finding the Darboux polynomials is always very costly computationally So, we 
have developed methods to extract those polynomials from the ODEs themselves. 

Some are very straightforward once we analyze the differential operator D = 
N d x + z Nd y + M d z , used on the procedure we have presented on [22] and, here, 
used as our example os methods benefitting from a quick finding of Darboux 
polynomials. Actually, as we have mentioned on the paper, this feature of the 
procedure developed on [22] is very useful. If N is a function of (x, y) it will be 
itself a Darboux polynomial and the same applies for any factor of iV that is a 
polynomial on (x,y). Analogous results apply for M (and its factors) that are 
polynomials on (z). Even this simple analysis allows that, in practice, one can 
use a Darboux type approach where a high degree Darboux polynomial is needed 
(say third degree and up) in a practical and feasible way. 

On the last sections of the paper, we have presented other examples of the 
idea of trying to extract the Darboux polynomials inspecting the "format" of the 
ODE. This time, a more elaborate analysis had to be performed. We have used 
the fact that the invariant is of the form: 



where A(x,y, z) ,T> (x,y, z) , B (x,y, z) and C(x,y,z) are all polynomials and, 
furthermore, T> (x, y,z) ,B (x, y, z) and C (x, y, z) are Darboux polynomials of 
the D-operator just mentioned. Of course, that is not the most general case (as 
previously stated). But, from our experience, this is not a very restricted case 
(in comparison to the stringent demand we have already imposed of considering 
only rational ODEs). From this we have the general expression for the ODE and, 
analyzing some classes of equations, were able to determine the needed Darboux 
polynomials for examples belonging to the class. 

With this, we hope to have establish that this kind of approach is very im- 
portant to the understanding of the structure of ODEs and to make the solving 
(reducing) of them, via Darboux methods, more practical in cases where high 
degree polynomials are present. 

Many roads are open to follow these ideas. We are currently pursuing some of 
them. For instance, we are working on a classificatory system for 20DEs based 
on equation (69) and the ideas here exposed. 
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